Asymptotic distribution of singular values 
of powers of random matrices 
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Abstract: Let i be a complex random variable such that Ex = 0, 
E\x\ 2 = 1, E\x\ 4 < oo. Let 

hi £ {1,2,...} be independet copies 
of x. Let X = (N~ 1 ^ 2 Xij), 1 < i, j < N be a random matrix. Writing 
X* for the adjoint matrix of X, consider the product x m X* m with some 
m £ {1, 2, . . . }. The matrix x m X* m is Hermitian positive semi-definite. 
Let Ai, Aa, . . . , Aat be eigenvalues of x m X* m (or squared singular values 
of the matrix X m ). In this paper we find the asymptotic distribution 
function 

G (m) {x)= lim EFj, m) (i) 

N— too 

of the empirical distribution function 

N 

F^(x) = N- 1 ^2l{X k <x}, 
k=l 

where I{^4} stands for the indicator function of event A. The moments 
of G (m) satisfy 

M' m) = f x p dG (m) (x) = 

Jr m-P- 

In Free Probability Theory Mp are known as Fuss-Catalan numbers. 
With m — 1 our result turns to a well known result of Marchenko-Pastur 
1967. 

Keywords and phrases: Random matrices, Fuss-Catalan numbers, 
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1. Introduction 

Let X = (A^ -1 / 2 ^-^), 1 < i, j < N be a random matrix. We assume that 
X{j = xVj are independent complex random variables such that 

B Xij = 0, E\ Xij \ 2 = l, Blxij] 4 < B (1.1) 
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with some B < oo independent of N. We assume additionally that 

L N (a)=N~ 2 E\ x ij\ Al {\ x ij\ > ctVN} -> as N -> oo (1.2) 

l<i,j<N 

for all a > 0. Note that Xj,- = and X = can depend on N, which 

is not reflected in our further notation. 

Writing X* for the adjoint matrix of X, consider the product 

■yy(m) _ ^mx*™ 1 

with some m E {1,2,...}. The matrix W"( m ) is Hermitian positive semi- 
definite. Let Ai, A2, . . . , \n be eigenvalues of W^" 1 ) (or squared singular val- 
ues of the matrix X m ). In this paper we find the asymptotic distribution 
function 

G^(x) = lim EF[ m) (x) 
of the empirical distribution function 

N 

F^\x) = N- 1 Y,H^k<x}, 
fe=l 

where I{^4} stands for the indicator function of event A. 

Theorem 1.1. Assume that (1.1) and (1.2) hold. Then the limit G^(x) = 
lirnzv-xx) E Fjf (x) exists. The function G^ m \x) is a distribution function 
and it has moments 

MW= [ x?dG^(x) = ^—( pm + P ). (1.3) 
Jr mp+l\ p J 

Corollary 1.2. Let Xij be independet copies of a random variable, say x, 
such that 

Ei = 0, E|x| 2 = l, E|x| 4 <oo. 

Let X = (N~ l l 2 Xij), 1 < i,j < N. Then the limit lirn2v-»-oo E Fpj (x) exists 
and it is equal to G^ m '(x). 

Gessel and Xin 2006 [4] showed that for any natural m the sequence 

. . . is a sequence of moments of some probability measure. Hence, 

G (m) 

is a probability distribution for any natural m. Since Mi m) < (? m with 
some c m < 00, by Carleman's Theorem in [3] the measure G^ is uniquely 
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determined by its moments. The support of the measure G^ is the interval 
[0,m- m (m + l) m+1 ]. 

With m = 1 Theorem 1.1 turns to a well known result of Marchenko- 
Pastur 1967 [6]. Namely, the asymptotic distribution of eigenvalues 

of the matrices XX* has the moments = ^pj(^f)- Note that in the 

case m = 1 our fourth moment assumption is stronger than assumptions in 
Theorem 2.5 and Theorem 2.8 in Bai 1999 [1]. The question of the weakest 
sufficient conditions in the case m > 1 remains an open problem. 

In Free Probability Theory Mp m ^ are known as Fuss-Catalan numbers. 
Combinatorial properties of this sequence have been studied by Nica and 
Speicher 2006 [8]. Mlotkowski 2009 [7] investigated a family of distribu- 
tions, say G^ m ' r \ with real m > and < r < m, such that G^ m,r ' has 
moments mp l p+r ( m P+ p+r ). It is easy to check that G^ = G (m \ Oravecz 
2001 [9] proved that powers of Voiculescu's circular element have distribu- 
tion G^ m \ This distribution belongs to the class of Free Bessel Laws (see 
Banica et al 2008 [2]). 

Let M. m (x) = Ylp^=o^p xP be the generating function of the sequence 

Mp m \ It satisfies the following functional equation (see equation (7.68) on 
p. 347 Graham et al 1988 [5]) 

M m (x) = l + xMZ +1 {x). (1.4) 

Equation (1.4) allows us to describe G^ in the framework of Free Proba- 
bility Theory. In Free Probability Theory the free multiplicative convolution 
£ M 1] is defined for any positive random variables £ and rj (see Nica and Spe- 
icher 2006 [8], p. 287). The S-transform is a homomorphism with respect to 
free multiplicative convolution, i.e. if £ and r] are free independent positive 
variables, then S^rj( z ) = S^(z)S r] (z). Recall that the 5-transform, say S(z), 
of a distribution \i is defined as follows. Let 

■ oo 

M p = I x p d/j(x), u(z) = y2M p z p . 
® p=i 

Then 

S(z) = —u-\z), (1.5) 

z 

where u~ l denotes the inverse function of u. 

Equation (1.4) allows to calculate the 5-transform, say S^ m \z), of G^ m \ 
and 

S( m \z) = - 1 - . (1.6) 
w (1 + z) m v ' 
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It means that the family G^ m ' has the following property: if a random vari- 
able £ has distribution then the r-th power of the S-transform of £ 
is equal to the S'-transform of multiplicative free power £ Kr . This property 
holds for this family of distributions only. 

To prove Theorem 1.1 we use truncation and the method of moments. 
Truncation means that we can replace (see Section 2.1 for details) X by the 
matrix X = (Xij) with truncated entries (here and below Xij = N~ 1 / 2 Xij 
denote entries of matrix X) 



= XijI{\Xtj\ < «at}, 



(1.7) 



where ajy is some sequence of positive numbers such that a at — > as N — > 
oo. Lemma 2.1 (see Section 2.1) reduces the proof of Theorem 1.1 to the 
proof of the next proposition. 



Proposition 1.3. Assume that ajy 
holds if 



and — > 0. Then Theorem 1.1 



X 



(N) 
ij 



< UN, 



max 

±<i,j<N 



EI 



(JV) 
ij 



< (3 N N- 3 / 2 , 



E 



X 



(JV) 
ij 



1/N 



<(3 N X 
(1.8) 



-3/2 



Let us explain our proof of Proposition 1.3. Denote by £ m (N) a random 



variable with distribution E F^ 1 ' . We show that the moments E £m{N) con- 
verge to Mp 171 ^ . In order to simplify the notation assume for a while that X^ 
are real random variables. Then one can represent ~E£ m (N) as 



2mp— 1 



3=0 



where the sum ^( 2mp ) i s taken over io, ..,ii mv G {1, .., N} such that i2mv = 
in. The notation Xfty means Xf ■ := Xj.j.,, in case of e(j) = + and 
Xj~£ := Xi j+1 i. in case of e(j) = — (see Section 2.2 for a precise definition 
of the spin variable e(j)). We investigate properties of paths (io, ~,i2 m p) by 
combinatorial methods. The moment E£^(iV) converges to the number of 
paths of a special type. Namely, one can describe such paths as follows: the 
cardinality of {io, .., i^mp) is equal to mp + 1 and each factor Xi^ +1 appears 
in the product rX/=o _1 X ti- +1 twice. I n Section 2.4 we count the number of 
these paths. 
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2. The proof of the main result. 
2.1. Truncation. 

Recalling that Xij = N~^/ 2 Xij, we can rewrite L^{a) as 

L N (a)= Yl E\ x ij\ 4 H\Xi,j\ > «}• 
i<i,j<N 

Since for all a > the ratio L^(a)/a i tends to 0, one can find a sequence 
ol jy J, such that L^{a^)j 'a A N — > and N 5 a^ oo for any 5 > as 

N — > oo. Let F^ 1 (t) denote the empirical spectral distribution function of 
the matrix X m X* m . 

Lemma 2.1. The limit behaviors ofEF^ n \t) and E F^\t) are the same, 
that is 

S up\BF^(t)-BFP(t)\^0. 

teR 

Proof. Since by definition l-Fjy 71 (i) — -F^(i)| 7^ only if there exist i, j € 
{1, . . . , iV} such that |Xy| > a^v, we have 

|EF^(t)-E4 m) (t)|< £ P(I*<*I>«n). (2- 1 ) 

l<i,j<N 

Estimating P(|Xjj I > a^r) < a^ 4 E |Xy-| 4 I{|Xy| > a^v} and using inequality 
(2.1) we obtain 

N 

S up\EF^\t)-EF^\t)\ < a^y j E\X l3 \H{\X l] \>a N } 

= a^L N (a N )^0. (2.2) 

□ 

Note, that the lower order moments of the truncated variables are asymp- 
totically equal to the moments of the original variables. Writing for a while 
X = X^ we have for k < 3 

\EX k - EX k \ < E\X\ k I{\X\ > a N }. (2.3) 
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The right hand side of (2.3) can be estimated as 

E |X| fc I{|X| > a N ] < c4~ 4 E \X\ 4 < /37viV~ 3/2 , (2.4) 

where (3 N = Ba^N" 1 ' 2 -> as N -> oo. 

Lemma 2.1 shows that the limit behaviors of F^(t) and F^(t) are the 
same. Thus we may replace X by X in the following arguments and assume 
that X is truncated, that is, that entries of X satisfy the assumption (1.8). 

2.2. Moments of the spectral distribution. 

We apply the method of moments. Recall that Ai, A2, • • • , Xn denote the 
eigenvalues of X m X* m . We can write 

TV 

E (iV) = iV^E X p j = iV -1 E Tr(X m X* m ) p . (2.5) 
j=l 

We assume that m and p are fixed and study the asymptotics of E£ m (iV) 
as N —7- 00. In order to simplify notation, hence forth we assume that Xy 
are real random variables. 

In the Hermitian case, the trace of X 2fc may be rewritten in terms of the 
entries of X via 

(2k) 2fc_1 

ETrX 2fe = ^ 'E [] (2.6) 

j=0 

where the sum * s taken over io, .., i s £ {1, .., N} such that i s = iq. 

In the non-Hermitian case E Tr(X m X* m ) p has a similar representation. 
An entry of X m X* m is given by 

P^- ]ik = 'y 1 X^Xi^ • • • Xi m _ 1 i m Xi m+1 i m • • • Xf : i 2m _ 1 (2.7) 

l<ij<N 

(2.8) 

We write Xf ■ := X; .« and XT.a := X; Then the right hand side 

ijtj+i 'j'j+i 'j+i'j 

of (2.8) takes the form 

2m- 1 

[x m n» fc = Yl II X $L> ( 2 - 9 ) 

l<ij<N j=0 

where iq = i, i2m = k, and the 'spin' variable e(j) takes values e(j) = + with 
j < m, and e{j) = - with j > m. Since (X m X* m )?> = X m X* m • • • X m X* m 
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(p times), one needs to change the order of indices in Xf^. +i if the spin 
e = — and 

£ (j) = j +,if -7 (mod 2m) E {0, ... m — 1}, 
I — , if j (mod 2m) G {m, . . . 2m — 1}. 

Using these notions (2.5) takes the form 

E£P(7V) = iV _1 E Tr(X m X* m )P 

2mp— 1 

= ^^v-e n *a>„- p.n) 

i=o 

A crucial notion in the proof is that of 'paths' of indices of the type (io, i\, . . . , i2mp~i, 
2.3. Description of paths. 

We consider a path i = (*0j ■ ■ ■ j*2mp-i) which corresponds to a product 

ngr 1 • Let ^ be a set ° f p airs {(i> i+i) £0) }i=c7 2 u{(2m P - 1, on, 

where (j,j + 1)+ := (j, j + 1), (j, j + 1)~ := (j + I J) and e(j) is given 
by (2.10). We call pairs + l) e ^ and (k, k + l) £(7c ) equivalent (denoted 

by (j,j + l) £ij) ~ (A;, fc + l) £ «) iff = We also call (j, j + l) e W 

an edge of the path i. We construct a directed graph Q\ as follows. A vertex 
V of Q\ is a subset of {0, 1, . . . , 2mp — 1} such that j £ V and k G V if and 
only if = There exists an edge (V,W) if and only if there exist I G V 
and r €U such that (Z, r) &V (note that | Z — t* | = 1). Denote by V the total 
number of vertices of the graph Q\ and by E its total number of edges. Since 
the graph Q is connected E > V — 1. It is clear that V is a cardinality of 
{io,«i, • • • ,«2m-i} and -E is a cardinality of a quotient set V / ~. Denote by 
k r (r = 1, . . . , E) the cardinality of each equivalence class in V . Note, that 
k\ + k 2 H Yk E = 2mp. 

Remark 2.2. Consider paths i = (io 3 ■ • • > hmp-i) and k = (&o, fci, . . . , &2m-i) 
such that Q\ = Q\a- It is clear that if identically distributed then 

2mp— 1 2mp— 1 

e n x S+i =e n x &+i- 

3=0 3=0 

We will show, that assuming our conditions the asymptotic products corre- 
sponding to equivalent paths are equal as well. 
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Definition 2.3. We define the contribution of a graph Q to 2.11 as 

2mp— 1 

Cont{G) = £ A^E II X $L 



i.Gi=g 



3=0 



Lemma 2.4. Using these notations we have that the contribution of the 
path Q is asymptotically given by 



(2.12) 



r=l 

when N tends to infinity. 
Proof. Since Xy are independent we have 

2mp— 1 _B 

e n v \i vx 



3=0 



r=l 



Furthermore, for any vertex V the number of possible values of corresponding 
indices (indices ij such that j E V) lies between N and N — 2mp ~ N. The 
lower bound iV — 2mp is due to the fact that indices corresponding to this 
vertex should not coincide with indices corresponding to other vertices and 
that there are at most Imp different indices. This yields the multiplicity N v . 
Together with the factor iV^ 1 this finally leads to the formula (2.12). □ 

Definition 2.5. We call a graph Q\ (m, p)-regular graph, if it has at least 
mp + 1 vertices and k r > 2 for all r G {1,2, . . . ,E}. The path i we call 
(m,p)-regular path. 

Lemma 2.6. Cont{Q\) does not converge to zero if and only if Q\ is the 
regular path. 



Proof. Since the variables Xij satisfy conditions (1.8), we have 



EI" 



(2.13) 



Of course, this estimaton holds for k = 1 too. At first we consider that one 
of k r is equal to 1 (without loss of generality k\ = 1). Then we have 



II E * 



i(r)j(r) 



r=l 



EX i(l)3(l) Tl EX 



i(r)j(r) 



r=2 



<p N N-^ 2 N- E+1 o^ Akr - 2) 



p N N-zl 2 N- E+1 a 2 ™ p - l - 2{E - 1] < N-z- 1 / 2 , 



(2.14) 
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and the contribution of such a graph is bounded by 

\Cont(i , i 2mp ^)\ < N v-i n -e-i/2 = N v-e-i n -V\ ( 2.i5) 

Note that V—E—l < since the graph Q is connected and hence N V ~ E ~ 1 N~ 1 / 2 
tends to 0. 

Furthermore, we consider the case V < mp + 1. Note that k r > 2 for any 
r and E < 2mp/2 = mp. Our truncation leads to 



r)j(r) 



r=l 



E ,E r (fcr-2) 



jV 



-E 2mp-2E 



N 



(2.16) 



Using inequality (2.16) to estimate the terms in (2.12), we obtain for such 
a product 



N 



v-i 



(r)j(r) 



r=l 



< N y - E ~ l a 2 ™ p - 2e . 



(2.17) 



Note that E > V — 1 and 2mp — 2E > 0. It follows that the right hand 
side of (2.17) does not converge to only if 2mp — 2E = and V — E — l = 0, 
i.e. V = mp + 1 and the graph Q is a regular graph. □ 

Furthermore, we obtain 

Lemma 2.7. A regular graph is a tree and it has exactly V = mp + 1 
vertices and exactly E = mp edges ( each representing an equivalence class 
of size k r = 2). 

Remark 2.8. Due to the fact that EX| ~ 1/N and by the remarks above 
we can write the contribution of a regular graph Q reg as 

Cont{g reg ) ~ 1. (2.18) 

We now show the connection between the moments of the spectral distri- 
bution "EiFm^ and the number of regular graphs. Indeed, ^ m (N) has distri- 
bution E Fm^ ■ Denote by T m ^ p the set of all possible graphs of view Q\ and 
by Tm,p the set of all (m,p)-regular graphs. Then 

EC(A0= E Cont(S)~ £ 1 = #T^. (2.19) 



We can reformulate 2.19 as 



Lemma 2.9. limjv->oo is equal to the number of (m,p) -regular graphs. 
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2-4- Counting of the number of regular graphs. 

Lemma 2.10. The number of all (m,p) -regular graphs is #T m e tP = M p m \ 
Proof. The numbers = ( mp p +p ) satisfy to the recurrence (see [5] ) : 

m— 1 

M ir ) ) = J2U M n ) ' M[ m) = l, (2.20) 

p-1 i=0 

where the sum Y^p-i ^ s taken over all po + p\ + • • • + p m = p — 1. We will 
show that there is one-to-one correspondence between collections of (m,pk)- 
regular graphs (Q m , po , G m , Pm ) ■ YZLoPi = P~ l and (ra, p)-regular graphs 
Gm,p- It follows that 

m m 

k% = # u ® t ^ = e n # r -i ( 2 - 21 ) 

p— 1 i:=0 p— 1 i:=0 

and the sequence #T m e iP satisfies to both the same reccurence and initial 
conditions as the sequence Fuss-Catalan numbers M p and, by this reason, 
these two sequences are equal. 

Proposition 2.11. The number if~T = 1 for all m. If Q\ is a (m, 1)- 

regular graph then indices ik and %i are equal iff (k + I) = 2m. 

Proof. By induction. Consider m = 1. In this case it is clear, that there 
is only one regular graph 0—^1 and Proposition 2.11 holds. Assume, that 
Proposition 2.11 holds for all m < mo- Consider the path i and a corre- 
sponding graph Q\ (see fig.l). This path has too + 1 distinct indices and 




Fig 1. The path i. Vertices, that correspond to equal indices, are connected via dotted lines. 

it has 2too at all. It follows that there exist at least 2 one-element ver- 
tices of Q\. Let these one-element vertices be {s} and {t}. Consider the 
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pair (i t -i , it) £ (t) . It must have an equal pair, but it is not equal to any 
other index. It means, that (if_i, itY^~^ = {h,H+i) £ ® ■ It follows, that 
e(t — 1) 7^ £ if)- There are exactly two possibilities for this: t = tuq or 
t = 0. Assume without loss of generality that s = and t = tuq. Therefore 
imo-l = imo+i (notice, that (mo — 1) + (mo + 1) = 2mo). Define (mo — 1, 1)- 
path j as follows: j k := i k if k G {0, . . . , m - 2}, j mQ _i := i mo _i = i mo+1 , 
jk '■= ik+2 if k G {mo,...,2(mo — 1) — 1}. (See fig. 2) The path j is the 




Fig 2. The path j is a (mo — 1, l)-path 



(mo — 1, l)-regular path. There is only one such path by inductive hypothesis 
and (i k = i{) <^ (j k = ji- 2 ) & (k + (l-2) = 2(m -l)) <^ (k + l = 2m ). □ 

Definition 2.12. Notice, that the vertex of a regular graph has two outgoing 
edges iff the corresponding index has the form iimk (because it should be 
(ij, ij+\) £ ^ = (ij,ij + \) and (ij-_i,ij) e w' _1 ) = and it happens if and 

only if j = 2mk). The distance between such vertex and vertex V is called a 
type of vertex V. The type of index ij is the type of a vertex V such that j G V. 
It is clear, that index ij has type j (mod 2m) if j (mod 2m) G {0, . . . , m— 1} 
or type —j (mod 2m) in the other case. There are m + 1 types of vertices. 
Note that only indices of the same type can be equal (this is proved in the 
case p = 1 in Proposition 2.11 and it will be proved for other cases below). 

Consider a collection of (m,pfc)-regular paths (io, ii, . . . , i m ) (such that 
J2T=oP k = P~^) an d collection of corresponding regular graphs (Qq, Q\, . . . , Q m ). 
Sum of path's lengths is 2m(p— 1). We indicate the recipe how to obtain the 
(m,p)-regular graph from these collections. We take an (m, l)-regular graph 
and attach to its vertices the graphs from the collection in the following 
way: the graph Qq is attached to vertex of type 0, the graph Q\ is attached 
to vertex of type 1, . . . , the graph Q m to the vertex of type m. For a more 
detailed argument we denote X^^=oP* by Pk (Pm = P — 1) and the indices of 
the k th path if, by ij . The resulting (m,p)-regular graph is denoted by Qy 
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Define the map A : A(io, ii, • • • , i m ) = j as follows 



.(o) . .(o) . .(o) 

JO ■= »o >J1 = H ^2mP -l : = l 2mp -V temP ■= JO] 

,-(!) „• _ ._.-(!) „• _ ._,-(!) 



j2mP +l 
j2mPi+2 
j2mP 2 +2 



"1 , • • • ,j2mP!-l := «2mpi-H J2mPi := «o ';J2mPi+l := j2mP +i; 
.(2) . .(2) . .(2) . .(2) 

h > • • • > J2mP 2 -l :— ^2mp2-l'-? 2mP 2 : ~ J ' J2mP 2 +l := h , 

j2mPi+2; 



_ -(m) ■ _ -(m) ■ _ -(m) 

j2mP m -i+m '■— lm , ■ ■ ■ , J2m(p-1)-1 : ~ l 2mp m —\ ' J2m(p— 1) :_ *0 ' 

J2mp— m+1 := j2mP m _i+m— li • • • > J2mp—k • = j2mPk+ki ■ ■ ■ i 32mp—l := j2mPi+l- 

(2.22) 

Let A is the corresponding map A(C/o,£/i) • • • ><5m) = Graphically the 
constructon (2.22) looks as follows: (fig. 3). 




Fig 3. The regular graph Gj, obtained from the collection of (m,pi) -regular graphs 
(Go, Gl, ■ ■ ■ , Gm) 



Example. For example, we consider for m = 2 the collection of (2,Pk)- 
regular graphs (£2,2 j $2,0, £2,1) (see fig. 4) and we obtain from it a (2,4)- 
regular graph ^2,4 (see fig.5). 

Proposition 2.13. Using the above construction we get an (m,p) -regular 
graph. 

Proof. Indeed, the graph Qj has exactly mp edges and k r = 2 for all r = 
1,2, ... , mp. Furthermore, there are exactly mp+1 vertices (there is no new- 
introduced vertex and there are exactly YaLo ( m Pi + 1) = m(p—l)+m+l = 
mp + 1 vertices of graphs Q^). □ 

Note, that the map A is the injection. 

Now we consider the arbitrary (m,p)-regular path i and try to construct 
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Fig 4. Graphs Q2.2 and C?2,i- Graph Q m ,o is empty 




Fig 5. The resulting graph C/2,4 



inverse map for A. Denote 
Jo := {j ■ ij = *o}; 

Jk ■= {j ■ 3 / 2mp - k, ij = i 2mp -k}, k G {1, . . . , m - 1}; ^ ^ 

J fc := max(J fc ), Jk := min(J fc ). 

We will prove that the sets Jk have some remarkable properties and after 
that it will be clear, how to obtain a collection of regular paths from one 
(big) regular path. 

Proposition 2.14. Jk is nonempty. 

Proof. Indeed, there is £ Jo and 2mp — m G J m . If Jk is void with 
k G {1, . . . , m — 1}, then the index i 2m p-k has no equal indices in the path 
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i. But in this case the pair (i^mp-k—i-i hmp-k)^ has no equivalent for the 
following reason. The index i2mv-k appears in (2mp — k, 2mp — k + 1) _ and 
in (2mp — k — 1, 2mp — k)~ only and they are not equivalent. But each edge 
in a regular path has equivalent one, a contradiction. Therefore the initial 
assumption that J/% is void must be false. □ 

Proposition 2.15. Jk (0 < k < m) are pairwise disjoint and if k < I then 
Jk < Ji (f° r o-H j ^ Jk an d f or M i £ J\ the inequality j < i holds). 

Proof. Indeed, if Jk f] Ji ^ then i 2 mp-k = ^imp-l- The edges of the path i 
have the same orientation on the section (2mp—k, 2mp—k—l)~ , . . . , (2mp— 
I + 1, 2mp — I) and therefore the graph G\ has a cycle. But a regular graph 
is a tree, a contradiction. Thus Jk f] Ji = 0- We prove the second part of 
Proposition 2.15 for the case I = fc+1 only (which is sufficient). Consider the 
edge (2mp — (k + 1), 2mp — k)~ . It must be equivalent to an edge (t, t + 1) + 
with some t G Jk and t+1 G Jk+i- If there exists s G Jk such that s > t then 
s > t + 1 (Jfc P| Jfc + i = 0). The edge (t,t + 1) is not equivalent to any edge 
in the section (t + l,t + 2), . . . , (s — 1, s), because it has only one equivalent 
edge (2mp — (k + 1), 2mp — k)~ . It follows that there are two different paths 
in the graph Gi which connect vertex U (such that t + 1 £ W) and vertex 
V (such that s G V and t G V), that is there is a cycle in the the graph Gi, 
and hence there is a contradiction. Therefore t = max Jk = Jk- Similarly, 
t + 1 = Jfc+i- It follows, that Jk + 1 = «/jfc-fi and for allj G </& and for all 
i G Jfc+i the inequality j < i holds. □ 

Proposition 2.16. For all k the difference (Jk — Jk) is divisible by 2m. 

Proof. Denote (J& — Jk) (mod 2m) by Notice, that (Jk — Jk) is the 
number of edges in the path's section ( Jfc, Jk + 1) , . . - , (Jfc — 1, Jk)- Notice, 
that the orientation of edges changes after every m steps. Edges of the 
form (ij-, ij k+1 ) + have the same orientation. It follows, that do < m — 1, 
d\ < m — 1, (do + 1 + cZi ) < m — 1 (mod 2m) (and so (do + 1 + d\) < m — 1, 
because < d + 1 + d\ < 2m - 1), . . . , < d + 1 + d\ + 1 + ... + d m - 2 + 
1 + d m —\ < m — 1 (similarly), i.e. < YlT=o dk + m — 1 < m — 1. Therefore, 
dk = for all k = 0, 1, . . . , m — 1. Consider all edges of the path i. There are 
m edges of the form (Jk,Jk+i), m edges of the form (2mp — k, 2mp — k + 1)~ 
with some k = 1,2, ... , m and all the remaing ones are in sections of the 
form (Jk, Jk + 1), • • • , (Jk ~ 1) ^fc)- There are 2mp edges in total. Therefore, 
SfcLo(^ — Jk)+m+m = 2mp and hence XlfcLo dk = ® (mod 2m). It follows, 
that d m = too. □ 
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Proposition 2.17. If Jk < t < Jk and J\ < s < J\, then i t / i s . In other 
words, sections of the path i of the form (Jk, J& + 1), . . . , (Jk — 1, Jk) with 
k = 0, 1, . . . , m are disjoint. 

Proof. Without loss of generality we consider I > k. Assume that it = i s . In 
this case the section (Jk, Jk + 1), • • • , (s — 1, s) contains the edge (Jk, Jk+i), 
and the section (t, t + 1), . . . , (Jk — 1, Jk) does not contain it or its equivalent 
(2mp — k—1, 2mp — k)~. Thus, there are two non-equal paths in the regular 
graph G\ which connected vertex hi (such that Jk E U) and vertex V (such 
that s£V and t E V), that is there is a cycle in the the graph G\. Therefore, 
the initial assumption must be false. □ 

Now we can describe the inverse map for A. Let pk '■= (Jk ~ Jk)/^ rn (Pk 
is a nonnegative integer by Proposition 2.16). Furthermore, we have for sum 
YlT=o Pk =p— 1 (see the proof of Proposition 2.16). Denote by j( fc ) the k-th 
resulting path (it has a length 2mpk and if pk = then is empty). Let 

3t '■= ijk+({t-k) mod 2m Pk )), t E {0, . . . ,2mp k — 1}, k E {0, ...,m}. (2.24) 

Now one obtains the collection (0^,2, {?2,o> ^2,1) (see fig. 4) from the graph 
£?2,4 (see fig.5) in the way described in (2.24). 

Proposition 2.18. The collection of paths (^°\ jW, . . . , j' m ^) (defined by 
(2.24)) is the collection of regular paths. 

Proof. In fact, the path is almost the same as the section (Jk,Jk + 
1), . . . , (Jk — 1, Jk) of the regular path i . This section contains 2mpk edges. 
Each of these edges has an equivalent one in the same section by Proposition 
2.17. Therefore this section contains exactly mpk + 1 distinct indices because 
of connectivity and non-cyclicity. Hence the path is a regular path. □ 

Thus, A is the bijection between Tm$ and |J Tm$ x Tm%\ x • • • x Tm$ m , 
where the union is taken over all po +pi + ■ ■ ■ +p m = p—1. Hence j^Tmf v = 
Mp™ 1 ^ and Lemma 2.10 is proved. □ 

Lemmas 2.9 and 2.10 show that the moments of the spectral distribution 
converge to Mp m \ Thus Theorem 1.1 is proved. 
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